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I.  INTRODUCTION 


GOALS  ASSOCIATED  WITH  MODEL  DEVELOPMENT 

Current  analysis  capabilities  for  determination  of  structural  response 
depend  greatly  on  the  ability  to  model  the  materials  involved.  The  adequacy 
of  future  system  design  and  analysis  therefore  depends  on  the  development  of 
advanced  material  models  for  reinforced  concrete.  The  goal  of  this  modeling 
effort  is  to  provide  a  more  complete  understanding  of  the  behavior  of  concrete 
and  to  provide  improved  capabilities  to  analytically  model  the  response  of 
reinforced  concrete  structures  to  applied  loads,  including  static  loads  and 
dynamic  blast  and  shock-induced  loads. 

OBJECTIVE  AND  SCOPE 

This  report  describes  a  continuing  evaluation  effort  on  the  hysteretic- 
fracturing  endochronic  concrete  model  developed  by  Bazant  and  Schieh  (Ref.  1) . 
This  model  is  intended  to  reproduce  all  significant  aspects  of  the  compressive 
behavior  of  concrete  except  strain  rate  effects.  Endochronic  models  assume 
that  inelastic  strain  accumulation  can  be  characterized  by  a  scalar  parameter 
called  intrinsic  time,  whose  increment  is  a  function  of  strain  increments.  In 
earlier  work  (Refs.  2  and  3)  it  was  shown  that  endochronic  theory  could  repro¬ 
duce  aspects  of  concrete  behavior  to  a  degree  never  before  experienced.  This 
model  is  a  refinement  of  the  earlier  work,  serving  to  improve  on  the  duplica¬ 
tion  of  observed  behavioral  features  while  removing  theoretical  inconsisten¬ 
cies.  Since  the  endochronic  model  is  assumed  to  be  one  for  the  compressive 
behavior  of  concrete  only,  it  was  meshed  with  a  directional  cracking  tension 
model  to  provide  a  complete  depiction  of  concrete  behavior. 

The  objective  of  tnis  study  is  an  evaluation  and  development  of  the 
hysteretic-fracturing  endochronic  model  as  a  viable  tool  for  the  analysis  of 
concrete  structures.  The  ability  of  the  model  to  synthesize  concrete  behavior 
for  a  variety  of  standard  and  nonstandard  load  paths  is  studied  through  com¬ 
parisons  made  with  the  triaxial  test  data  on  plain  concrete  recently  generated 
at  New  Mexico  State  University.  This  comparison  necessitated  the  investiga¬ 
tion  of  various  parameter  adjustment  techniques  for  fitting  the  model  to  con¬ 
cretes  outside  the  range  for  which  the  model  was  calibrated.  The  original 
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algorithm,  which  was  devised  to  operate  in  a  two-dimensional  finite  element 
structural  analysis  program,  is  generalized  to  three  dimensions  in  this  effort 
in  anticipation  of  its  insertion  into  a  more  genercl  code. 

Explosively  loaded  systems  respond  with  high  strain  rates  and  it  is 
known  that  concrete  is  a  strain-rate-sensitive  material.  Current  concrete 
models  do  not  attempt  to  account  for  this  effect.  Although  strain  rate 
effects  are  not  included  in  the  endochronic  concrete  model,  the  intrinsic  time 
term  provides  an  attractive  avenue  for  the  incorporation  of  such  effects. 
Possible  approaches  to  a  strain  rate  dependency  for  the  model  are  explored. 

One  aspect  of  the  endochronic  concrete  algorithm  which  is  often  cited 
as  a  major  disadvantage  is  the  large  number  of  parameters  which  must  be  stored 
for  each  element  for  each  increment  of  a  finite  element  analysis.  It  is 
important  to  note  that  constitutive  relations  which  attempt  to  provide  a  real¬ 
istic  depiction  of  hysteresis  require  a  large  number  of  storage  parameters,  so 
comparisons  between  the  endochronic  model  and  other  models  which  omit  hystere¬ 
sis  are  inappropriate.  It  is  apparent,  however,  that  the  number  of  parameters 
associated  with  the  model  must  be  reduced  if  the  model  is  to  be  suitable  for 
practical  application.  In  this  study,  modification  of  the  hysteresis  portion 
of  the  algorithm  is  considered  as  a  possible  approach  for  reduction  of  the 
nunber  of  required  element  storage  parameters. 
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II.  REVIEW  OF  THE  ENDOCHRONIC  CONCRETE  MODEL 

REVIEW  OF  THE  8AZANT-SCHIEH  ENDOCHRONIC  MODEL 

The  endochronic  model  utilized  in  this  analysis  consists  of  two  parts: 
a  directional  cracking  model  used  to  describe  tensile  behavior,  and  a  modified 
version  of  the  Bazant-Schieh  hysteretic-fracturi ng  endochronic  model  used  to 
describe  compressive  behavior.  The  Bazant-Schieh  model  was  intended  to  serve 
as  a  refinement  of  the  endochronic  concrete  model  developed  in  an  earlier 
work.  The  primary  improvements  to  the  previous  theory  involve  (1)  separate 
considerations  of  inelastic  strains  related  to  plastic  slip  of  the  micro- 
cracked  material  and  inelastic  strains  related  to  strength  and  rigidity 
degradation  due  to  microcracking,  (2)  a  more  detailed  modeling  of  inelastic 
volumetric  behavior,  and  (3)  institution  of  a  load-unload  criterion  along  with 
jump-kinematic  hardening  to  better  match  hysteresis  and  ensure  energy  dissipa¬ 
tion  under  cyclic  loading. 

The  intrinsic  time  parameters  used  in  this  model  are  defined  primarily 
in  terms  of  the  deviatoric  strain  path  length,  whose  increment  is 


where 

r.  is  the  deviatoric  strain  path  length 

de,  is  the  deviator  of  the  strain  tensor 
km 

A  separate  intrinsic  time  parameter  is  used  for  plastic  and  fracturing  devi¬ 
atoric  inelastic  strains,  so  that 

deii  =  — 11  ^  (1) 

J  2G 

and 


★ 

e.  .d< 
ij 


(2) 
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where 


de^j,  de^j  are,  respectively,  the  plastic  and  fracturing  inelastic 
deviatoric  strain  increments 

G  is  the  elastic  shear  modulus 

★ 

S.  .  =  S-  -  u.  . 
ij  ij  ij 

S..  is  the  deviatoric  stress 
ij 
★ 

e  ■  .  =  e  -  a! 
ij  ij  '  J 

e..  is  the  deviatoric  strain 
1 J 

a-j,  a|j  are  jump-kinematic  hardening  parameters 

d4  =  cF idi  (plastic  intrinsic  time  increment)  ( 

dk  =  cF^dt;  (fracture  intrinsic  time  increment) 

c  is  a  coefficient  defining  the  loading  condition 

Fl  and  F2  are  invariant  functions  of  stress,  strain,  and  5 


The  particular  forms  of  Fi  and  F2  are  determined  empirically.  Expressi 
similar  to  Equations  1  and  2  are  developed  for  the  dilatational  and  shei 
compaction  aspects  of  inelastic  volumetric  behavior.  However,  inelastic  volu¬ 
metric  behavior  due  to  pore  crushing  is  related  to  the  volumetric  strain  path 
length  rather  than  deviatoric  strain  path  length.  The  volumetric  strain  path 
length  is  defined  as 


The  earlier  endochronic  model  did  not  accurately  reproduce  the  hysteretic 
behavior  which  occurred  during  unload-reload  cycles  and,  in  fact,  violated 
Drucker's  stability  postulate.  Rectifying  these  problems  required  the  insti¬ 
tution  of  a  loading-unloading  criterion  with  approx  -iate  adjustments  for  these 
situations.  Separate  criteria  were  suggested  for  deviatoric  and  volumetric 
behavior  based,  respectively,  on  deviatoric  and  volumetric  work.  Adjustments 
to  unload-reload  cycles  are  then  partially  made  by  selecting  appropriate  val¬ 
ues  of  the  constant  parameter,  c,  in  Equations  3  and  4,  along  with  a  similar 
parameter  involved  in  the  pore  crushing  contribution  to  the  inelastic  volumet¬ 
ric  strains.  In  addition,  jump-kinematic  hardening  is  included  to  complete 
the  unload-reload  model. 


MODIFICATIONS  TO  THE  BAZANT-SCHIEH  MODEL 


The  parameters  defining  functions  Fj  and  F2  in  Equations  3  and  4  and 
similar  functions  for  inelastic  volumetric  behavior  were  determined  by  fitting 
data  to  a  form  arrived  at  from  a  combination  of  physical  reasoning  and  experi¬ 
ence.  The  resulting  expressions  are  relatively  complicated,  and  are  not 
unique--! .e. ,  variations  of  the  form  and  parameters  in  the  functions  could 
produce  equally  satisfactory  fits  to  the  experimental  data.  In  this  study  an 
algorithm  for  the  endochronic  theory  was  developed  and  encoded  to  facilitate 
the  investigation  of  the  model.  A  slight  variation  of  the  form  of  function  F2 
was  necessary  to  produce  a  satisfactory  fit  to  experimental  data.  The  result¬ 
ing  model  produced  results  which  differed  slightly  from  those  of  Reference  1. 

Some  questions  arose  concerning  the  jump-kinematic  hardening  associated 
with  the  model.  The  inelastic  deviatoric  strains  are  calculated  by  the 
equations 


(5) 


(6) 


as  described  previously.  At  the  onset  of  unloading  and  reloading,  the  a  terms 
and  the  c  term  change  values,  based  on  a  deviatoric  work  criterion 


dW  >  0,  W  =  WQ  :  virgin  loading  ;  c  =  1 

dW  <  0  :  unloading  ;  c  =  cy  =  0.6  ;  a—  =  S^-,  a)j  =  e^  (unloading) 

dW  >  0,  W  <  WQ  :  reloading  ;  c  =  cr  =  0.8  ;  ai  •  =  s^, 
a|  ■  =  e ;  •;  ( reloading) 


A  similar  formulation  is  used  for  volumetric  unloading  and  reloading. 

Although  this  approach  works  quite  well  when  the  material  is  completely 
unloaded,  it  appears  to  break  down  under  a  partial  unloading.  Since  the  a  and 
a'  values  in  Equations  5  and  6  are  defined  at  the  start  of  unloading  and  at 
the  start  of  reloading,  but  not  at  renewed  virgin  loading,  artificial  harden¬ 
ing  is  produced  whenever  the  material  is  not  subject  to  a  significant  unload¬ 
ing.  Accounting  for  this  artificial  hardening  was  apparently  a  consideration 
in  the  application  of  jump-kinematic  hardening  to  the  plastic  fracturing  con¬ 
crete  theory  of  Bazant  and  Kim  (Ref.  3).  A  remedy  that  appears  to  work  for 

this  situation  is  to  reset  the  a.,  and  a’.,  terms  to  zero  when  the  virgin  load- 

ij  ij  y 

ing  path  is  renewed.  A  similar  adjustment  is  made  for  volumetric  behavior. 

The  resulting  depiction  of  the  concrete  behavior  is  consistent  with  that 
observed  exper imental ly . 

REVIEW  OF  THE  DIRECTIONAL  CRACKING  MODEL 

The  endochronic  concrete  model  was  designed  to  reproduce  characteristic 
compressive  behavior  of  concrete.  Tensile  (cracking)  behavior  must  be  modeled 
separately  if  the  material  model  is  to  be  considered  for  general  finite  ele¬ 
ment  applications.  A  proper  cracking  model  works  with  principal  stresses,  can 
accommodate  compressive  loading  of  the  crack,  produces  gradual  cracking,  and 
allows  for  a  partial  resistance  to  shear  remaining  on  the  crack  plane. 

A  suitable  model  was  obtained  by  modifying  the  cracking  approach  used 
for  an  earlier  endochronic  formulation.  This  model  is  built  around  the  con¬ 
cept  that  cracking  produces  a  free  surface  so  that,  for  instance,  a  plane 
stress  problan  is  reduced  to  a  uniaxial  stress  formulation  through  the  forma¬ 
tion  of  a  crack  plane  normal  to  the  plane  of  stress.  Strains  perceived  by  the 
material  in  the  direction  normal  to  the  crack  are  not  the  strains  defined  by 
nodal  point  motion  in  a  finite  element  program,  but  instead  are  dependent  on 
the  material  behavior  in  the  direction  parallel  to  the  crack.  This  distinc¬ 
tion  is  crucial  for  continued  proper  compressive  modeling  using  the 


endochronic  model.  Crack  closure  is  then  defined  as  occurring  when  the  strain 
perceived  by  the  material  becomes  as  large  (compressi vely)  as  strains  defined 
by  nodel  point  motion. 

Instantaneous  crack  formation  for  an  element  in  a  finite  element  pro¬ 
gram  leads  to  a  phenomenon  known  as  crack  shock,  which  is  a  misrepresentation 
of  the  material  behavior  due  to  the  inability  of  the  finite  element  program  to 
properly  ha-mile  discontinuous  behavior.  Crack  shock  can  be  avoided  by  allow¬ 
ing  the  crack  plane  to  develop  into  a  free  surface  gradually,  rather  than 
instantaneously.  This  is  done  by  gradually  dropping  the  tensile  stress  down 
to  zero.  The  transition  condition  is  similar  to  the  free  surface  condition, 
except  the  stress  on  the  crack  surface  has  some  specified  value  other  than 
zero. 

The  principal  modifications  made  to  the  earlier  cracking  model  involved 
introduction  of  a  principal  stress  formulation  and  inclusion  of  a  partial 
resistance  to  shear  on  the  crack  plane.  The  initial  crack  check  is  made  by 
considering  stresses  in  the  principal  directions.  It  is  known  that  some 
reduced  shear  resistance  exists  on  a  rough  crack  surface.  However,  retention 
of  shear  resistance  on  crack  faces,  along  with  the  assumption  that  cracks 
occur  on  principal  planes,  implies  a  second  crack  need  not  be  perpendicular  to 
the  first  crack.  A  model  which  accounts  for  this  behavior  would  be  quite 
complicated.  For  the  sake  of  simplicity,  it  was  assumed  that  a  second  crack 
be  normal  to  the  first  crack.  This  tensile  model  produced  satisfactory 
results  for  test  problems  and  direct  applications. 


III.  PATH-DEPENDENT  CHARACTER  I  ST  ICS  OF  THE  ENDOCHRON1C  CONCRETE  MODEL 


INTRODUCTION 

Path  dependence  is  produced  in  the  endochronic  model  through  the  defi¬ 
nition  of  the  intrinsic  time  parameters  in  terms  of  the  strain  path  length. 
Whether  or  not  the  intrinsic  time  term  accurately  models  damage  accumulation 
during  the  loading  process  must  be  determined  through  a  consideration  of 
experimental  data. 

In  an  earlier  pnase  of  this  effort  (Ref.  4),  a  number  of  stress  paths 
were  considered  to  determine  what  stress  path  effects  would  be  predicted  by 
the  model  and  to  determine  whether  some  failure  surface  could  be  associated 
with  the  failure  conditions  predicted  by  the  thoery.  This  study  revealed  that 
(1)  the  model  implies  a  failure  surface  in  the  deviatoric  plane  which  is  con¬ 
sistent  in  form  with  that  observed  experimentally,  and  (2)  the  nature  of  the 
endochronic  model  is  to  produce  failure  dependence  on  the  stress  path,  other 
than  what  would  be  predicted  by  a  failure  surface. 

The  monotonical ly  increasing  aspect  of  the  intrinsic  time  term  in  endo¬ 
chronic  theory  would  cause  it  to  produce  inelastic  strain  increments  for  all 
deviatoric  strain  increments,  including  loading  tangent  to  what  could  be  con¬ 
sidered  the  failure  surface.  The  resulting  predicted  behavior  would  be  some¬ 
what  softer  than  observed  concrete  behavior  and  considerably  softer  than  the 
elastic  behavior  predicted  by  a  classical  plasticity  approach  to  inelastic 
behavior  (Ref.  5).  This  indicates  that  the  endochronic  theory  is  more  conser¬ 
vative  than  classical  plasticity  in  its  depiction  of  loading  to  the  side. 

COMPARISON  WITH  EXPERIMENTAL  DATA 

Calibration  of  mode1--In  this  portion  of  the  report,  comparisons  are 
made  of  predictions  by  the  endochronic  concrete  model  for  failure  conditions 
and  stress  path  dependence  with  tri axial  data  developed  recently  (Ref.  6). 

The  test  program  at  New  Mexico  State  University  (NMSU)  was  intended  to  provide 
information  on  hardening,  limit  surfaces,  path  dependency,  and  the  nature  of 
inelastic  deformation  in  general.  A  set  of  standard  compressive  stress  paths 
was  run  using  cubical  specimens,  followed  by  a  series  of  nonstandard  stress 
paths  intended  to  indicate  stress  path  dependency. 


The  endochronic  concrete  model  can  be  adapted  to  many  concretes  simply 
by  applying  the  appropriate  uniaxial  cylinder  compressive  strength  (f£). 
However,  this  method  was  unsatisfactory  for  the  NMSU  concrete,  which  had  a 
relatively  low  strength  (25.5  MPa)  in  comparison  to  the  concretes  considered 
for  determination  of  the  internal  parameters  for  the  endochronic  model.  Con¬ 
sequently,  an  approach  recommended  by  the  authors  of  the  theory  was  used  to 
calibrate  the  model.  After  some  preliminary  runs  to  study  the  effect  of 
selected  internal  parameters  on  the  calculated  response,  adjustments  were  made 
to  produce  as  good  a  fit  as  possible  to  the  uniaxial,  biaxial,  and  triaxial 
test  data.  The  adjusted  parameters  were  then  used  in  the  model  for  the  dup¬ 
lication  of  the  other  stress  paths  considered  in  the  test  series.  Matchups  to 
uniaxial  and  biaxial  behavior  are  shown  in  Figures  1  and  2  respectively. 

Several  comments  must  be  made  concerning  the  adjustment  of  parameters 
for  materials  other  than  those  to  which  the  model  had  originally  been  cali¬ 
brated.  Uniaxial  compressive  strength  alone  does  not  provide  a  complete  des¬ 
cription  of  a  concrete;  other  factors  such  as  water  content  can  cause  signifi¬ 
cant  variations  in  behavior.  The  concrete  used  in  the  NMSU  test  series  dis¬ 
played  some  features  that  are  dissimilar  to  those  of  the  concretes  used  to 
define  the  original  parameters.  For  instance,  shear  enhanced  compaction  was 
noticeably  reduced,  while  dilatation  was  observed  to  a  higher  degree.  Also, 
the  uniaxial  strengths  of  cubical  specimens  are  commonly  lower  than  those  of 
cylindrical  specimens.  Because  parameter  definition  in  the  model  is  based 
primarily  on  specimen  behavior  that  is  similar  to  the  behavior  of  cylinders, 
the  uniaxial  strength  of  cylindrical  specimens  was  used  in  the  calibration. 
Finally,  because  the  parameters  were  adjusted  by  trial  and  error,  the  number 
of  parameters  varied  was  kept  to  a  minimum.  The  model  may  be  more  sensitive 
to  some  parameters  than  to  others;  adjustments  of  this  type  may  therefore  have 
undesirable  rami f icat i ons  for  load  paths  other  than  those  matched.  More 
recently  a  curve-fitter  program  has  been  used  for  these  adjustments.  This 
program  allows  a  larqer  number  of  parameters  to  be  considered. 

Comparison  of  stress-str ai n  rel ations--Predictions  based  on  the  endo¬ 
chronic  model  were  compared  to  the  experimental  data  from  two  perspectives: 

(1)  stress-strain  plots  were  used  to  illustrate  the  nature  of  the  inelastic 
deformation  produced  by  the  model,  and  (2)  the  predicted  and  measured  failure 
surfaces  were  evaluated.  Because  features  evident  in  some  stress-strain  plots 
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are  duplicated  in  others,  not  all  the  plots  from  the  stress  paths  considered 
are  displayed. 

Calculated  and  experimental  data  are  compared  in  Figure  3  for  a  loading 
in  which  pure  shear  is  superimposed  on  a  hydrostatic  loading.  Agreement  is 
generally  good,  although  the  model  somewhat  underpredicts  inelastic  behavior 
for  the  unloaded  and  neutrally  loaded  directions. 

A  biaxial  stress  state  with  a  stress  increment  in  the  second  principal 
direction  that  is  one-half  that  in  the  first  principal  direction  is  shown  in 
Figure  4.  The  loading  in  Figure  4  is  repeated  in  the  data  shown  in  Figure  5, 
with  the  addition  that  the  loading  in  the  third  principal  direction  is  equal 
to  one-tenth  the  loading  in  the  first  principal  direction.  The  strength 
enhancement  due  to  this  relatively  small  confinement  effect  is  visible  in  both 
the  calculated  and  the  experimental  data.  The  higher  stiffness  observed  in 
the  calculated  data  may  be  related  to  the  use  of  third  stress  invariant  terms 
within  the  model  for  adjusting  to  triaxial  loadings.  It  is  possible  that 
these  terms  overadjust  for  cases  close  to  a  biaxial  condition  because  similar 
difficulties  were  not  apparent  for  more  balanced  triaxial  loadings.  The 
parameter-adjustment  technique  used  in  this  part  of  the  study  would  exaggerate 
this  effect. 

Three  load  paths  were  considered  which  were  variations  of  the  biaxial 
loading  state  shown  in  Figure  2.  Figure  6  displays  the  behavior  resulting 
from  a  loading  path  consisting  of  a  hydrostatic  load  to  0.1  f£  followed  by 
an  equal  biaxial  load.  The  strength  enhancement  produced  by  this  relatively 
small  amount  of  confinement  can  be  seen  by  comparison  with  Figure  2  and  is 
suitably  depicted  by  the  model. 

A  cycle  is  added  to  the  hydrostatic  loading  for  the  behavior  depicted 
in  Figure  7.  A  sliqht  hardening  is  evident  in  both  the  experimental  and  the 
calculated  data  as  a  result  of  the  cyclic  loadinq. 

The  behavior  displayed  in  Figure  8  is  the  result  of  a  loading  path  con¬ 
sisting  of  a  hydrostatic  load  to  0.1  f^  followed  by  a  uniaxial  load-unload 


Strain 


Figure  4.  Stress  versus  strain  for  the  proportional  stress  path  02 


Figure  5.  Stress  versus  strain  for  the  proportional  stress  path 
a2  =  0.5oi,  o 3  ■  Q.loi. 


Figure  6.  Stress  versus  strain  for  an  equal  biaxial  loading  superimposed 
on  a  hydrostatic  state. 
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Figure  7.  Stress  versus  strain  for  a  loading  consisting  of  a  hydrostatic 
state  followed  by  a  partial  hydrostatic  unload,  and  then  an 
equal  biaxial  loading. 
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Figure  8.  Stress  versus  strain  for  a  loading  consisting  of  a  hydrostatic 
state  followed  by  a  single  uniaxial  cycle,  and  then  an  equal 
biaxial  loading. 
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cycle  and  then  an  equal  biaxial  loading  to  failure.  The  general  match  to  the 
experimental  data  is  good  but  the  model  does  not  produce  the  amount  of  hyster¬ 
esis  found  in  the  experimental  data.  It  is  convenient  to  compare  the  behavior 
shown  in  Figure  8  with  that  shown  in  Figure  6,  where  the  results  of  a  similar 
loading,  but  with  the  cycle  omitted,  are  displayed.  An  apparent  increase  in 
strength  is  seen  in  the  experimental  data  for  the  load  that  includes  cycling. 
This  behavior  is  exaggerated  in  the  model.  A  possible  explanation  relates  to 
the  approach  used  to  produce  hysteresis  in  the  model.  When  deviatoric  unload¬ 
ing  begins,  the  value  of  the  deviatoric  strain  path  length  is  stored.  The 
deviatoric  strain  path  length  is  reset  to  this  value  when  the  unload-reload 
cycle  ends,  as  defined  by  the  restoration  of  the  deviatoric  strain  energy  lost 
during  unloading.  (This  adjustment  is  necessary  to  prevent  an  exaggeration  of 
hardening  in  the  model  due  to  accumulation  of  the  deviatoric  strain  path 
length  during  cyclic  loading.)  It  is  suspected  that  the  change  in  the  nature 
of  the  loading  causes  the  deviatoric  strain  path  length  to  be  reset  prema¬ 
turely,  which  produces  higher  effective  values  for  the  path  length  than  is 
appropriate  and  leads  to  an  overly  stiff  behavior  for  the  model.  This  hypoth¬ 
esis  obtains  credence  from  the  behavior  shown  in  Figure  7,  where  a  hydrostatic 
cycle  is  substituted  for  the  uniaxial  cycle.  No  excessive  strength  enhance¬ 
ment  occurs  for  the  model  in  this  case  because  the  hydrostatic  cycle  does  not 
produce  a  deviatoric  strain  path  length. 

Similar  comparisons  were  made  for  combinations  of  hydrostatic  and  uni¬ 
axial  loads.  In  Figure  9,  behavior  is  depicted  for  a  loading  path  consisting 
of  a  hydrostatic  load  to  0.1  r,  followed  by  a  uniaxial  loading  to  failure. 
Comparison  with  Figure  1  demonstrates  the  strength  enhancement  due  to 
confinement  present  in  both  the  experimental  and  calculated  data.  The  model 
is  seen  to  significantly  underpredict  the  strength  of  the  material  under  these 
conditions.  This  behavior  is  apparently  related  to  the  sensitivity  of  the 
model  to  the  invariant  terms  discussed  previously,  since  the  phenomenon  was 
not  observed  in  the  calibration  triaxial  test  where  the  initial  hydrostatic 
load  was  about  0.4  f\  Figure  10  shows  the  behavior  for  a  similar  loading, 
but  an  equal  biaxial  loading  cycle  is  interjected  between  the  hydrostatic 
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Figure  9.  Stress  versus  strain  for  a  loading  consisting  of  a  hydrostatic 
state  followed  by  a  uniaxial  loading  to  failure. 
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Figure  10.  Stress  versus  strain  for  a  loading  consisting  of  a  hydrostatic 
state  followed  by  an  equal  biaxial  cycle,  and  then  a  uniaxial 
loading  to  failure. 


and  uniaxial  loadings.  This  cycle  is  seen  to  provide  a  strength  enhancement 
for  both  the  experimental  and  calculated  data,  but  the  calculated  data  under¬ 
predicts  the  experimental  data  to  approximately  the  same  degree  as  seen 
earlier. 

Results  from  a  loading  path  consisting  of  a  hydrostatic  load-partial 
hydrostatic  unload,  and  then  a  pure  shear  load  to  failure  are  depicted  in  Fig¬ 
ure  11.  Comparison  with  Figure  3  which  depicts  the  same  loading  without  the 
hydrostatic  cycle  reinforces  the  conclusion  previously  drawn  that  hydrostatic 
cycling  in  this  load  range  has  little  influence  on  material  behavior. 

The  last  of  the  more  complicated  compressive  loading  paths  to  be  dis¬ 
cussed  provide  some  insight  into  path-dependent  characteristics  of  the  mate¬ 
rial  behavior.  Figure  12  displays  results  for  a  loading  path  consisting  of  a 
hydrostatic  load  to  0.1  f^,  followed  by  a  uniaxial  loading  to  0.85  f ' , 
followed  by  a  shear  load  until  the  stress  condition  in  the  direction  of  the 
uniaxial  load  is  switched  with  that  of  one  of  the  other  loaded  directions, 
followed  by  a  uniaxial  load  to  failure.  This  rather  intricate  loading  path 
involves  a  120-deg  excursion  in  the  deviatoric  plane,  as  shown  by  path  A  in 
Figure  13.  Such  a  loading  path  would  not  produce  significant  inelastic  behav¬ 
ior  for  a  material  which  followed  a  classical  plasticity  law.  Figure  9,  which 
depicts  a  similar  loading  path,  but  with  no  excursion,  provides  a  convenient 
datum.  It  can  be  seen  that  while  the  excursion  results  in  a  slight  reduction 
in  the  strength  for  the  experimental  data,  a  significant  hardening  occurs  for 
the  calculated  data.  Such  behavior  is  difficult  to  interpret,  but  appears  to 
be  related  to  the  discontinuous  behavior  which  can  result  from  the  hysteresis 
portion  of  the  model.  The  deviatoric  strain  path  length  upon  which  inelastic 
behavior  is  based  is  redefined  if  an  unloading  cycle  occurs,  so  loading  paths 
(such  as  this  one)  wnich  are  very  close  to  a  neutral  loading  can  produce 
ambiguous  results. 

Figure  14  depicts  a  loading  path  consisting  of  an  equal  biaxial  load  to 

0.1  f'  followed  by  a  uniaxial  load  to  0.75  f',  and  then  a  uniaxial  load  in  the 

c  c 

second  direction  to  0.75  f^,  followed  by  a  biaxial  load  to  failure.  This 
loading  path  examines  the  effect  of  a  loading  which  roughly  parallels  the 
failure  line  in  a  yrJj  -  I,  invariant  space.  A  classical  plasticity 
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Figure  12.  Stress  versus  strain  for  a  loading  consisting  of  a  hydrostatic 
state  followed  by  a  uniaxial  loading,  and  then  a  change  of 
principal  axes  followed  by  a  uniaxial  load  to  failure. 
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Figure  13.  Loading  path  of  Figure  12  in  deviatoric  stress  space 
(tension  positive). 
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approach  would  be  expected  to  display  little  inelastic  behavior.  A  convenient 
loading  path  for  comparison  is  shown  in  Figure  2  which  displays  an  equal  biax¬ 
ial  loading.  The  two  paths  are  depicted  in  deviatoric  space  as  paths  B  and  C, 
respectively,  in  Figure  9.  A  slight  reduction  in  strength  is  seen  for  the 
experimental  data,  but  a  strength  increase  is  predicted  by  the  theory.  It  is 
expected  that  the  increased  deviatoric  strain  path  length  artificially  hardens 
the  calculated  behavior  depiction  for  the  more  complicated  path.  Apparently 
the  presence  of  inelastic  behavior  for  such  paths  is  countered  in  the  endo- 
chronic  theory  by  the  hardening  effect  produced  by  the  longer  deviatoric 
strain  path. 

In  Figure  15,  the  following  loading  path  is  depicted:  a  hydrostatic 
loading  to  2  f^;  then  a  shear  cycle  of  Aoj  =  ao2  =  f^;  then  a  hydrostatic 
unload  to  0.1  f'  followed  by  an  equal  biaxial  loading  to  failure.  The  three 
specimens  used  for  this  loading  path  showed  widely  varying  results,  possibly 
because  of  the  difficulty  involved  in  following  a  prescribed  loading  path  at 
high  stresses.  Comparison  with  Figure  6,  which  shows  the  results  of  a  hydro¬ 
static  loading  to  0.1  f^  followed  by  an  equal  biaxial  loading,  indicates  a 
slight  strength  gain  for  the  experimental  data.  The  model  produces  an  exag¬ 
gerated  strength  gain,  probably  as  a  result  of  the  accumulation  of  the  devi¬ 
atoric  strain  path  length.  Otherwise,  the  model  produces  results  within  the 
spread  of  the  experimental  data.  The  kinks  which  occur  in  the  calculated 
results  are  due  to  the  beginning  or  ending  of  volumetric  or  deviatoric  inelas¬ 
tic  behavior. 

COMPARISON  WITH  STRESS  PATHS  INVOLVING  TENSION 

Loading  paths  involving  tension  were  also  included  in  the  New  Mexico 
State  University  study.  Although  the  results  of  these  tests  are  useful  for 
the  evaluation  of  the  total  concrete  model,  the  tensile  cracking  model  domi¬ 
nates  behavior  rather  than  the  endochronic  (compressive)  model.  It  is  impor¬ 
tant  to  note  that  the  stress  levels  and  strain  levels  reached  in  these  tests 
are  quite  low.  This  causes  a  magnification  of  the  influence  of  procedure- 
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Stress  versus  strain  for  a  loading  consisting  of  a  hydrostatic 
state  followed  by  a  shear  cycle,  and  then  a  partial  hydrostatic 
unload,  and  then  an  equal  biaxial  loading  to  failure. 


related  phenomena,  such  as  specimen  rotation,  which  adversely  effect  the  accu¬ 
racy  of  the  data.  A’so  note  that  the  assumption  of  isotropic  behavior  in 
concrete  is  an  approximation.  Anisotropy  had  an  extreme  influence  on  several 
tests,  particularly  the  equal  biaxial  tensile  loading  path.  Consequently,  the 
only  loading  paths  discussed  are  those  which  provide  clear  insight  into  the 
behavior  of  the  model. 

A  uniaxial  tension  loading  path  is  depicted  in  Figure  16.  It  is  impor¬ 
tant  to  emphasize  that  the  stress  and  strain  levels  are  an  order  of  magnitude 
lower  than  those  commonly  occurring  for  compressive  loads.  The  theoretical 
concrete  cracking  strength,  defined  using  a  standard  formula  (Ref.  7),  gave  a 
good  representation  of  the  strength  of  the  material.  The  model  predicted  a 
stiffen  response  than  was  observed  in  the  experiment.  Anisotropic  behavior 
appears  in  the  experimental  data  for  the  lateral  direction  but  is  not  produced 
in  the  isotropic  moael.  Nonlinear  behavior  is  seen  in  the  experimental  data 
but  is  not  obvious  in  the  calculated  data. 

A  loading  path  consisting  of  a  uniaxial  load  to  0.5  f£  followed  by 
a  shear  load  until  failure  is  described  in  Figure  17.  The  experimental  data 
indicate  failure  at  a  stress  slightly  lower  than  the  f£  value.  Failure 
for  the  model  is  defined  as  attainment  of  a  stress  equal  to  f{.  The  cal¬ 
culated  and  experimental  data  match  well.  The  kink  seen  in  plots  of  both  sets 
of  data  occurs  when  the  loading  changes  from  uniaxial  stress  to  pure  shear. 

Study  of  comparisons  of  predicted  versus  experimental  data  for  other 
tensile  loading  paths  reinforce  the  following  observations:  (1)  Path  depend¬ 
ency  for  the  failure  stress  is  apparent  in  the  experimental  data  but  is  not 
reproducible  in  the  model;  (2)  The  calculated  data  appear  to  match  the  experi¬ 
mental  data  quite  well,  when  test-procedure-related  effects  are  discounted; 

(3)  The  model  had  a  tendency  to  give  a  slightly  stiffer  representation  of  the 
material  behavior,  as  was  also  seen  in  compression  at  low  stresses. 

SUMMARY 

The  failure  points  for  those  compressive  loading  paths  in  which  failure 
occurred  at  a  lode  angle  of  either  0  deg  or  60  deg  are  plotted  in  /J2  -  I, 
invariant  stress  space  in  Figure  18.  A  linear  regression  fit  matched  the 
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Figure  16.  Stress  versus  strain  for  a  uniaxial  tension  loading  path. 
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Figure  17.  Stress  versus  strain  for  a  shear  loading  superimposed  on  a 
uniaxial  loading. 


experimental  and  calculated  data  well  for  the  load  ranges  considered.  No 
attempt  was  made  at  a  similar  comparison  for  the  tensile  loading  paths  because 
the  maximum  principal  stress  cracking  criterion  does  not  result  in  a  surface 
in  this  space. 

When  any  concrete  model  is  evaluated  by  comparing  it  with  experimental 
data,  one  must  consider  the  high  degree  of  variation  in  concrete  behavior 
caused  by  different  aggregates,  water  contents,  and  design  strengths.  The 
method  used  in  this  study  for  parameter  adjustment  may  also  influence  the 
calculated  response.  However,  several  general  observations  can  be  drawn  from 
the  comparisons  of  the  experimental  and  calculated  data: 

1.  The  model  provides  a  good  qualitative  match  to  the  behavioral  fea¬ 
tures  observed  in  the  material. 

2.  The  quantitative  match  of  the  model  to  the  experimental  data  varied 
but  in  general  was  quite  good. 

3.  Quantitative  matching  of  the  response  features  suffers  in  certain 
cases  when  unloading-reloading  cycles  occur. 

4.  Any  adjustment  of  the  model  to  materials  outside  the  range  consid¬ 
ered  in  fitting  the  initial  parameters  must  be  made  with  caution  so  that  the 
reproduction  of  general  stress  paths  is  not  influenced  unduly. 

5.  Path  dependence  in  the  material  is  generally  exaggerated  in  the 
model . 
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IV.  GENERALIZATION  OF  THE  ENDOCHRONIC  ALGORITHM  TO  THREE  DIMENSIONS 

INTRODUCTION 

A  three-dimensional  model  was  developed  using  the  basic  ideas  previously 
incorporated  for  use  in  the  endochronic  constitutive  model  as  implemented  in 
the  two-dimension  finite  element  code  SAMS0N2  (Ref.  8).  The  resulting  model 
can  be  used  in  two-dimensional  codes  by  setting  the  appropriate  snear  compon¬ 
ents  of  the  stress  and  strain  tensor  to  zero. 

Generalizing  the  two  dimensional  endochronic  model  to  three  dimensions 
involves  three  steps:  (1)  upgrading  appropriate  arrays  and  loops  to  the  more 
general  case;  (2)  choosing  and  implementing  a  three-dimensional  eigenvalue- 
eigenvector  solution  technique;  and  (3)  implementing  a  general  second  order 
tensor  axis  transformation  scheme.  Only  the  last  two  steps  require 
elaboration. 

DETERMINATION  OF  PRINCIPAL  STRESSES  AND  DIRECTIONS 

Given  a  state  of  stress  at  a  point  defined  by  the  stress  tensor  o  with 

respect  to  a  basis  (e  ,  e  ,  e  ]  there  exist  three  mutually  orthogonal  planes 

x  y  z 

on  which  the  shear  stresses  vanish.  These  are  the  principal  planes,  and  the 
normal  stresses  acting  on  these  planes  are  the  principal  stresses.  It  is  the 
principal  stresses  on  which  the  cracking  criteria  are  based.  The  problem  to 
be  solved  is  the  eigenvalue-eigenvector  problem  expressed  in  direct  tensor 
notation  as 

(2  -  *  l)  *  X  =  0 

where  X  is  an  eigenvalue  or  principal  stress  and 
principal  direction  associated  with  X.  ^  is  the 
system  of  linear  homogeneous  equations  to  have  a 
minant  of  the  coefficients  must  vanish. 

det  -  »5iJ]  -  0 


(7) 

X  is  the  eignevector  or 
identity  tensor.  For  a 
nonzero  solution,  the  deter- 

(8) 


where 


6ij  = 


i  =  j:  5=1 


i  *  j:  6=0, 


The  result  of  performing  the  operation  defined  in  Equation  8  results  in 
a  cubic  equation  in  X  which  may  be  solved  using  any  of  a  number  of  standard 
techniques.  The  solution  method  chosen  here  is  explicit  and  requires  no  iter¬ 
ation.  It  follows  the  method  described  in  Reference  9. 


Expressing  the  determinant  in  terms  of  the  deviatoric  stress  invariants 


gives 


where 


(X-  )3  -  J'  X*  -  J'  =  0 
2  3 


x1  =  any  one  of  the  principal  deviator  stresses 

,.e.,  r  =  »,  -iokk 


with  repeated  index  implying  summation: 


J'  =  second  invariant  of  the  deviatoric  stress  tensor 
2 

=  l  [^n  -  °22)2  +  (°22  "  a33)2  +  (°33  "  an)2] 


+  °122  +°232  +  °3l‘ 


J'  =  third  invariant  of  the  deviatoric  stress  tensor 
3  =  det  (ai j • ) 


where 


°i j  =  °ij  "  3  6ij  °kk 
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Using  the  substitution 


X1  =  2  j;)1/2  cos  a 


in  Equation  9  gives 

2  /—  [4  COS3  a  -  3  cos  a]  =  J'  (10) 

\3  V  3 

Note  that 

4  COS3  a  -  3  COS  a  =  cos  3a 
Equation  10  becomes 


J'  /•>  \  3/2 

COS  3a  =  _i  (L.)  (11) 

2  l  J1  J 

Let  aj  satisfy  0  <  3aj  <  v  in  Equation  11  and  note  that  3aj,  3a!  +  2",  and 
3«!  -  2n  all  have  the  sane  cosine;  then  the  desired  expression  for  the  princi¬ 
pal  deviatoric  and  total  stresses  are  obtained  with 


where 


2n 

a3  =  <*l  - 

3 
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The  principal  directions  are  computed  by  solving  for  ^  in  Equation  7  for  each 

computed  A  and  then  normalizing  to  unit  length.  Note  that  the  occurrence  of 

principal  stress  of  multiplicity  greater  than  one  is  possible.  For  three 

equal  principal  stresses  (hydrostatic  stress  state)  the  principal  directions 

are  set  to  the  global  (e  ,  e  ,  e  )  directions.  For  two  equal  principal 

~x  ~y 

stresses  the  directions  associated  with  the  two  different  stresses  are  com¬ 
puted  and  the  third  is  computed  using  the  cross-product  rule. 

TRANSFORMATION  TO  PRINCIPAL  OR  CRACKED  BASIS 

Once  the  principal  stresses  and  directions  have  been  computed  the 
components  of  the  stress  and  strain  tensors  are  transformed  from  the  global 
basis  to  the  principal  or  cracked  basis  system.  This  is  performed  by  the 
transformation 

afo  =  aAi  aBj  °ij 

where  are  the  components  of  the  tensor  referred  to  the  new  basis,  and  a^1 
is  the  transformation  matrix  whose  components  are 


where  e?  is  the  component  of  unit  vector  i.  If  these  are  the  principal 
basis  vectors  the  transformat  ion  results  in  a  diagonal  matrix  whose  entries 
are  the  principal  stresses.  Checking  the  principal  stress  against  the  tensile 
limit  determines  whether  a  crack  is  formed.  Once  a  crack  is  noted  the  direc¬ 
tion  associated  with  the  crack  is  recorded.  When  only  one  crack  exists  it  is 
necessary  to  determine  the  other  two  directions  orthogonal  to  the  crack  direc¬ 
tion  which  reduce  the  shear  stresses  parallel  to  the  crack  plane  to  zero.  In 
other  words,  one  seeks  the  directions  orthoqonal  to  the  crack  plane  on  which 
the  maximum  normal  stresses  exist.  When  two  cracks  exist,  the  third  direction 
is  specified  due  to  the  assumption  that  cracks  form  orthogonal  to  each  other. 
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Once  the  basis  vectors  (whether  principal  or  crack  direction)  are 
known,  the  transformation  is  performed.  Thus  the  constitutive  model  works  in 
the  principal  basis  when  no  cracks  are  present  and  in  the  crack  direction 
basis  if  there  are  cracks  (even  when  the  cracks  have  closed).  On  exit  from 
the  constitutive  law  the  inverse  transformation  is  performed  so  that  global 
stresses  and  strains  are  available  for  printout  and  other  required 
computations. 

Figure  19  is  a  flow  chart  of  steps  involved  in  the  directional  cracking 
analys i s . 
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Figure  19.  Flow  chart  of  the  directional  cracking  model. 
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V.  HIGH  STRAIN  RATE  EFFECTS  IN  THE  ENOOCHRONIC  CONCRETE  MODEL 

INTRODUCTION 

The  impulse  loads  considered  for  the  analysis  of  protective  structures 
are  of  such  a  magnitude  that  significant  strain  rate  effects  are  expected  in 
the  material  response.  This  effect  influences  the  structural  response  through 
a  strength  enhancement  effect  at  high  strain  rates  and  through  resulting 
effects  on  stress  distribution  and  redistribution  in  indeterminate  structures. 
Existing  concrete  models  either  ignore  these  effects  or  treat  them  in  an  arti¬ 
ficial  fashion.  In  this  section  a  modified  version  of  the  endochronic  con¬ 
crete  model  which  deals  with  strain  rate  effects  is  described.  This  f^proach 
is  simple  and  flexible,  but  is  in  a  form  which  appears  compatible  with  the 
physical  phenomena  involved. 

REVIEW  OF  CONCRETE  STRAIN  RATE  EFFECTS 

The  current  state  of  knowledge  in  the  area  of  strain  rate  effects  on 
the  compressive  behavior  of  concrete  is  summarized  in  recent  reports  by  Bazant 
and  Oh  (Ref.  10)  and  Suaris  and  Shah  (Ref.  11).  A  study  of  these  summaries 
and  other  specific  reports  on  experimental  studies  (Refs.  12  and  13)  reveals 
two  interesting  points  about  the  current  perception  on  strain  rate  effects  on 
concrete.  First,  most  of  the  existing  experimental  data  are  concerned  with 
testing  under  nominally  uniaxial  stress  conditions.  Second,  strains  are  often 
measured  in  an  indirect  fashion  which  creates  some  question  about  any  conclu¬ 
sion  about  performance  at  a  given  strain  rate.  Despite  the  nebulous  nature  of 
the  existing  data,  some  conclusions  are  currently  accepted  with  regard  to  high 
strain  rate  effects.  It  appears  that,  with  an  increase  in  strain  rate,  the 
extent  of  internal  mi crocr ack i ng  decreases,  resulting  in  less  nonlinearity  in 
the  stress-strain  relations  at  high  strain  rates.  The  result  is  an  increase 
in  strength,  secant  modulus,  and  strain  at  failure  for  the  concrete.  Although 
these  conclusions  are  primarily  based  on  uniaxial  stress  conditions,  similar 
results  were  seen  in  the  limited  amount  of  multiaxial  data  available. 
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ENOOCHRONIC  APPROACH  TO  HIGH  STRAIN  RATE  EFFECTS  IN  CONCRETE 


Although  strain  rate  effects  were  considered  in  the  earlier  development 
of  the  endochronic  concrete  theory  described  in  Reference  2,  the  concern  was 
creep  rather  than  high  strain  rate  effects.  An  endochronic  approach  has  been 
applied  to  metals  subjected  to  a  constant  strain  rate  by  Lin  and  Wu  (Ref.  14), 
who  later  developed  a  formulation  applicable  to  variable  strain  rates  in 
metals  (Ref.  15).  This  work  serves  to  indicate  the  potential  of  the  endo¬ 
chronic  model  for  representing  strain  rate  effects.  Specifically,  the  intrin¬ 
sic  time  term  provides  an  attractive  avenue  for  inclusion  of  strain  rate 
effects . 

Consider  a  heuristic  example:  A  representative  intrinsic  time  incre¬ 
ment  term  could  be  of  the  form 

dc  =  F  j (o,e)(dekmdekm) 1/2 
or,  alternately 

dA  =  Fj  (o  ,e )  (  dt 

where 

(  *  )  indicates  differentiation  with  respect  to  time. 

Neither  of  the  above  expressions  has  a  rate  dependency.  However,  a  variation 
of  the  form 


d;  =  F1(c,e)(ekmekm  +  C)i/2  dt 

would  have  a  rate  dependency.  The  C  term  could  represent  a  constant  or  a 
function.  The  intrinsic  time  increment  term  could  be  rewritten  as 


dc  =  Fj(o,e 


1  + 


1/2 


'kmekm . 


■  ^ekm^ekm 


1/2 


whicn  is  the  same  as  the  first  expression  considered,  except  for  the  brack- 
etted  term,  which  introduces  strain  rate  effects.  For  this  simple  model  it  is 
seen  that  low  strain  rates  produce  larger  intrinsic  time  increment  terms  than 
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high  strain  rates.  Since  the  intrinsic  time  term  serves  as  a  measure  of  dam¬ 
age  (microcracki ng)  in  the  concrete,  an  adjustment  of  this  form  can  be  related 
to  the  decrease  in  internal  microcracking  hypothesized  from  experimental  data. 
The  reduction  of  the  intrinsic  time  increment  in  the  model  reduces  the  inelas¬ 
tic  portion  of  the  material  response,  resulting  in  a  strengthening  of  the 
materi  a  1 . 

Since  different  intrinsic  time  terms  are  used  in  the  endochronic  con- 
ciete  model  for  plastic  and  fracturing  deviatoric  behavior  and  for  some  ine¬ 
lastic  volumetric  behavior,  different  strain  rate  considerations  could  be  made 
for  these  different  phenomena.  The  need  for  such  considerations  was  investi¬ 
gated  using  data  from  uniaxial  testing  with  the  corresponding  average  strain 
rate  taken  as  constant.  Applying  the  high  strain  rate  modification  only  to 
the  intrinsic  time  term  related  to  plastic  deviatoric  behavior  produced  the 
most  satisfactory  results.  Consequently,  the  formulation  described  in  this 
section  is  constrained  in  that  fashion.  Subsequent  availability  of  more 
detailed  experimental  data  may  require  an  adjustment  to  this  approach. 

The  simple  formulation  described  above  is  unsuitable  for  general  use  if 
C  is  a  constant,  but  serves  to  indicate  the  ability  of  the  theory  to  accommo¬ 
date  strain  rate  effects.  A  more  appropriate  expression  could  be  determined 
by  choosing  the  C-term  to  be  a  function.  A  convenient  alternative  is  to 
replace  the  entire  bracketted  term  with  a  term  consistent  with  observed  high 
strain  rate  behavior.  A  suitable  starting  point  is  a  modification  of  an 
expression  used  for  metals  for  a  similar  purpose  in  Reference  15: 

dc  =  k(dekmdekm)1/2 

where 


ka,  kb,  e0  are  parameters  to  be  determined 


The  influence  of  different  k  values  on  the  response  of  the  endochronic 
model  is  described  in  Figure  20  for  uniaxial  stress  loadings.  Increased 
strength,  secant  modulus,  and  strain  at  failure  are  all  reflected  in  the 
behavior  of  the  model. 

In  evaluating  the  validity  of  the  expression  given  above  for  k,  as  in 
determining  appropriate  parameters  for  the  expression,  one  is  limited  by  the 
sparsity  of  experimental  data.  Operating,  once  again,  within  the  constraint 
of  uniaxial  stress  data  and  utilizing  the  average  strain  rate  as  if  it  were 
constant,  one  can  define  appropriate  values  for  the  k  ,  kb,  and  eo  parameters 
by  using  summary  data  such  as  that  presented  in  Figure  21  (Ref.  11).  The 
resulting  parameters  are  not  unique  but  appear  to  produce  reasonable  results. 
They  are 

ka  =  0.09 
kb  =  0.0015 
e0  =  0.5  E-5 


SUMMARY 

Concrete  strain  rate  effects  are  reviewed  in  this  section  and  a  modifi¬ 
cation  to  the  endochronic  concrete  model  to  reflect  these  effects  is  proposed. 
The  modification  ties  strain  rate  effects  into  the  portion  of  the  model  which 
deals  with  the  physical  phenomena  considered  to  cause  such  effects. 

The  resulting  model  is  not  appreciably  more  complicated  than  the  time- 
independent  endochronic  model  and  is  amenable  to  change  if  such  change  appears 
warranted  by  advances  in  knowledge  of  strain  rate  effects.  A  detailed  check¬ 
out  of  the  model  through  comparison  with  existing  multidimensional  test  data 
is  desirable. 
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Figure  21.  Effect  of  strain  rate  on  concrete  strength  (from  Ref.  11) 


VI.  REDUCTION  OF  STORAGE  REQUIREMENTS 


INTRODUCTION 

The  present  form  of  the  endochronic  algorithm  requires  that  27  parame¬ 
ters  be  stored  per  element  for  each  load/time  step.  Storage  requirements  of 
this  magnitude  are  commonly  required  for  models  which  include  hysteresis 
effects.  Indeed,  elimination  of  the  hysteresis-related  stored  parameters 
would  leave  only  10  stored  parameters  for  the  combined  endochronic  and  crack¬ 
ing  model.  Nevertheless,  the  capability  of  the  model  as  an  effective  analysis 
tool  is  inhibited  by  excessive  storage  requirements.  Although  inclusion  of 
hysceresis  effects  makes  the  algorithm  more  rigorous,  such  rigor  may  be  neces¬ 
sary  only  when  the  model  is  used  as  a  research  tool  or  for  problems  involving 
significant  cycling.  Consequently,  simplification  of  the  approach  used  for 
cyclic  loading  is  an  appropriate  consideration  for  attaining  a  reduction  in 
the  storage  requirements  for  the  model. 

ALTERNATE  SCHEMES  FOR  HYSTERESIS 

The  approach  currently  used  for  hysteresis  requires  6  storage  elements 
to  define  the  loading  states  of  unloading/reloading/virgin  loading,  2  storage 
elements  for  adjustments  to  the  equations  based  on  the  loading  state,  and  9 
storage  elements  to  describe  jump-kinematic  hardening.  Three  alternate 
schemes  were  considered  for  simplifying  the  model.  The  first  was  complete 
elimination  of  specific  hysteresis  modeling  within  the  algorithm.  This 
approach  eliminates  the  17  storage  parameters  related  to  hysteresis,  but  would 
produce  unadjusted  endochronic  unloading  and  reloading,  a  condition  which 
implied  the  creation  of  energy  in  earlier  endochronic  models.  Excessive  accu¬ 
mulation  of  the  deviatoric  strain  path  length  would  also  be  expected  to  occur 
under  cyclic  loading,  resulting  in  an  artificial  hardening  of  the  material. 

The  second  scheme  considered  was  intended  to  be  a  compromise  between 
the  total  elimination  of  hysteresis  and  the  existing  model.  In  this  version 
the  adjustment  which  prevented  excessive  accumulation  of  the  deviatoric  strain 
path  length  was  retained  from  the  original  model.  It  is  necessary  to  differ¬ 
entiate  between  virgin/renewed-virqin  loading  and  unloading/reloading  for  this 
approach.  The  number  of  storage  parameters  related  to  hysteresis  is  reduced 
to  6. 


The  third  approach  which  was  studied  used  an  elastic  description  for 
unloading  and  reloading.  A  means  of  differentiating  between  virgin/renewed 
virgin  loading  and  unloading/reloading  is  again  necessary.  The  endochronic 
portion  of  the  algorithm  is  bypassed  during  unloading  and  reloading,  so  no 
incremental  inelastic  behavior  is  produced.  A  modest  increase  in  efficiency 
would  be  expected  as  a  consequence.  This  scheme  results  in  a  reduction  to  4 
of  the  number  of  storage  parameters  related  to  hysteresis. 

EVALUATION  OF  ALTERNATE  SCHEMES 

Evaluation  of  alternatives  to  the  complete  hysteresis  model  presently 
in  the  model  was  done  by  considering  only  uniaxial  loading  paths,  since  most 
of  the  available  data  on  hysteretic  behavior  fall  into  this  category.  Behav¬ 
ior  for  multidimensional  loadings  can  be  extrapolated  for  each  scheme,  to  some 
degree,  from  the  uniaxial  paths. 

One  loading  path  considered  was  a  cyclic  loading  consisting  of  complete 
unloading/complete  reloading  cycles.  This  is  a  duplication  of  a  path  for 
which  experimental  data  are  available  (Ref.  16).  The  matchup  to  these  data 
for  the  complete  hysteresis  model  is  shown  in  Figure  22.  The  model  produces 
hysteresis,  follows  the  general  envelope  of  the  experimental  data,  and 
reflects  the  degradation  of  the  elastic  modulus. 

Results  for  the  alternate  scheme  in  which  hysteresis  modeling  was  elim¬ 
inated  are  compared  to  the  experimental  results  in  Figure  23.  Hysteresis  is 
eliminated;  in  fact,  Drucker's  postulate  is  violated.  Degradation  of  the 
elastic  modulus  is  seen  to  occur.  The  model  appears  to  produce  a  reasonable 
depiction  of  the  envelope  of  the  experimental  data,  although  it  overshoots  the 
envelope  in  the  high  strain  range. 

S'.ie  insight  into  the  behavior  described  above  can  be  found  by  studying 
Figure  24  which  compares  the  same  experimental  results  with  results  from  the 
second  approximate  scheme.  This  approach  included  the  adjustment  which  pre¬ 
vented  excessive  accumulation  of  the  deviatoric  strain  path  length.  Since  it 
has  been  shown  that  the  excessive  accumulation  of  the  deviatoric  strain  path 
length  results  in  hardening,  it  was  expected  that  elimination  of  this  accumu¬ 
lation  would  result  in  a  somewhat  softer  response;  however,  the  extreme  soft¬ 
ening  requires  an  explanation. 


Figure  22.  Stress  versus  strain  for  uniaxial  cyclic  compression  using  the 
original  hydrostatic  endochronic  algorithm. 


23.  Stress  versus  strain  for  uniaxial  cyclic  compress 
endochronic  algorithm  without  hysteretic  consider 
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Figure  24.  Stress  versus  strain  for  uniaxial  cyclic  compression  using  the 
modified  hysteretic  endochronic  algorithm. 


The  problem  considered  was  cycled  by  loading  to  a  given  strain,  and  then 
unloading  to  a  zero  stress.  In  the  original  model  the  stiffness  associated 
with  reloading  is  lost  when  virgin  loading  is  renewed  for  two  reasons:  (1) 
the  center  for  jump-kinematic  hardening  is  redefined  (not  applicable  in  this 
case),  and  (2)  adjustments  to  the  formulation  are  made,  in  particular,  the 
deviatoric  strain  path  length  is  reset.  Renewed  virgin  loading  is  defined  to 
occur  when  the  deviatoric  strain  energy  lost  in  the  unloading  process  is 
restored.  When  the  hysteresis  model  is  completely  eliminated,  one  would 
expect  softer  behavior,  since  the  strain  at  which  the  next  unload  cycle  should 
occur  would  be  reached  at  a  lower  stress  level  without  the  influence  of  hys¬ 
teresis.  The  reason  the  first  scheme,  which  excluded  special  consideration 
for  unloading,  did  not  appear  softer  is  because  excessive  accumulation  of  the 
deviatoric  strain  path  length  produced  a  compensating  hardening  effect.  The 
modified  version  which  did  not  allow  such  excessive  accumulation  produced  a 
much  softer  and  unacceptable  response.  Since  the  interaction  of  these  two 
effects  in  the  first  scheme  cannot  be  expected  to  be  so  accommodating  for 
general  loading  conditions,  one  must  suspect  that  the  first  scheme  may  also  be 
undesirable  as  a  simplifying  approximation.  It  should  be  emphasized,  however, 
that  either  scheme  could  be  suitable  for  a  loading  which  is  known  to  be  free 
of  significant  unloading. 

The  third  approach  considered,  which  utilized  an  elastic  scheme  for 
unloading  and  reloading,  is  compared  with  the  experimental  data  in  Figure  25. 
Degradation  of  the  elastic  modulus  can  be  seen,  but  hysteresis  is  not  repre¬ 
sented.  The  envelope  of  the  exper oriental  data  is  violated  because  renewed 
virgin  loading  is  defined  at  a  higher  stress  than  for  the  other  schemes  due  to 
the  elastic  behavior  during  the  cycle. 

Generally,  a  simplified  version  of  the  hysteresis  modeling  would  be 
utilized  for  problems  where  hysteresis  would  not  be  expected  to  have  a  signif¬ 
icant  influence.  However,  there  is  evidence  that  the  hysteresis  model  is  a 
sensitive  aspect  of  the  overall  endochronic  formulation.  For  instance,  the 
modification  described  in  Section  I  of  this  report  resulted  from  a  slight 
fault  in  the  jump-kinematic  hardening  model  which  produced  an  inordinant 
amount  of  hardening.  This  behavior  is  displayed  in  Figure  26,  where  a  virgin 
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Figure  25.  Stress  versus  strain  for  uniaxial  cyclic  compression  using  the 
endochronic  algorithm  with  elastic  cycling. 


Cyclic  load 


hardening  representation. 


loading  curve  is  compared  to  a  similar  loading  which  included  a  small  unload- 
reload  cycle.  A  similar  comparison  is  made  in  Figure  27,  where  the  modified 
version  of  the  model  is  used.  Factors  such  as  the  virgin  loading  definition 
still  have  some  influence  on  the  results,  but  the  representation  is  accept¬ 
able.  It  is  important  to  note  that  behavior  such  as  this  is  not  uncommon  for 
concrete,  where  a  loading  redistribution  due  to  cracking  or  strain  softening 
could  result  in  a  slight  unloading  at  a  point  in  the  material.  The  first  and 
third  schemes  were  studied  under  this  type  of  load.  In  Figure  28  the  virgin 
loading  curve  is  superposed  in  the  results  of  the  first  scheme,  which  elimi¬ 
nated  hysteresis  considerations,  for  the  load  which  included  the  small  cycle. 
The  energy  stored  in  this  model  is  exaggerated  compared  to  the  virgin  loading. 
This  excess  is  due  to  accumulation  of  the  deviatoric  strain  path  length  during 
the  cycle.  A  much  more  acceptable  comparison  is  seen  in  Figure  29,  where  the 
scheme  with  the  elastic  unloading  is  considered. 

Another  concern  with  simplified  approaches  to  hysteresis  is  whether  the 
model  is  consistent  with  the  virgin  loading  curve  as  an  envelope  during  multi¬ 
ple  cycles.  Figure  30  compares  this  envelope  with  the  existing  model.  Some 
violation  of  the  envelope  is  seen  tc  occur.  The  loading  is  repeated  in  Fig¬ 
ure  31  for  the  alternate  scheme  which  ignores  hysteresis.  The  effect  of  an 
accumulation  of  the  intrinsic  time  term  is  obvious  in  the  extreme  violation  of 
the  envelope.  Similar  comparison'  are  made  in  Figure  32  for  the  scheme  which 
unloaded/reloaded  elastically.  This  comparison  is  severely  affected  by  the 
loading  technigue.  The  loadings  were  applied  until  a  given  strain  was 
reached,  while  the  unloading  was  stopped  at  a  given  stress.  Thus,  the  results 
were  affected  by  the  lack  of  hysteresis  representation  in  the  model.  A  simi¬ 
lar  misrepresentation  could  occur  for  a  general  loading.  Note  that  similar 
problems  would  exist  for  any  model  that  unloads  and  reloads  elastically. 

SUMMARY 

Any  of  the  alternate  versions  of  the  endochronic  model  will  give  rea¬ 
sonable  results  if  no  hysteresis  is  present.  The  scheme  in  which  hysteresis 
considerations  were  eliminated  tends  to  produce  exaggerated  hardness  in  the 
presence  of  cycling.  The  balance  between  softness  due  to  endochronic  unload¬ 
ing  and  hardness  due  to  accumulation  of  deviatoric  strain  path  length,  which 
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Figure  27.  Stress  versus  strain  for  uniaxial  compression  with  a  single  partial 
unloading  cycle  using  the  modified  jump-kinematic  hardening 
representation. 


Cyclic  load 


Figure  28.  Stress  versus  strain  for  uniaxial  compression  with  a  single  part 
unloading  cycle  for  the  endochronic  algorithm  without  hysteretic 
considerations. 
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Figure  29.  Stress  versus  strain  for  uniaxial  compression  with  a  single  partial 
unloading  cycle  for  the  endochronic  algorithm  with  elastic  cycling. 
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Stress  versus  strain  for  uniaxial  compression  with  multiple  cycli 
using  the  endochronic  algorithm  without  hysteretic  considerations 


Cyclic  load 


32.  Stress  versus  strain  for  uniaxial  compression  with  mi 
using  the  endochronlc  algorithm  with  elastic  cycling 


was  seen  for  t tie  uniaxial  cyclic  loading,  should  not  be  expected  to  exist  for 
general  loads.  The  approach  which  adjusted  for  the  accumulation  of  deviatoric 
strain  path  length  responded  much  too  softly  in  the  presence  of  unloading. 

The  scheme  which  defined  unloading  and  reloading  elastically  appeared  to  give 
the  most  reasonable  results. 


VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


CONCLUSIONS 

The  endochronic  concrete  algorithm  as  considered  in  this  study  contin¬ 
ues  to  gain  credibility  as  a  viable  concrete  model.  It  demonstrates  a  capa¬ 
bility  of  reflecting  all  significant  behavioral  characteristics  of  concrete 
under  static  loading  and  has  shown  the  flexibility  required  for  the  incorpora¬ 
tion  of  dynamic  effects.  Conclusions  derived  from  this  effort  are 

1.  The  model  provides  a  good  qualitative  match  to  the  behavioral  fea¬ 
tures  observed  in  the  material  for  a  variety  of  loading  paths.  This  perfor¬ 
mance  lends  credence  to  the  ability  of  the  model  to  handle  general  loading 
paths. 

2.  The  model  fares  worst  under  circumstances  for  which  the  intrinsic 
time  term  accumulates  excessively.  This  was  particularly  evident  for  loadings 
which  paralleled  the  failure  surface. 

3.  Adjustment  of  the  model  to  concretes  outside  the  range  considered 

in  calibrating  the  original  model  should  be  done  with  caution  to  avoid  unduly 
influencing  the  reproduction  of  general  stress  paths. 

4.  The  most  promising  approach  to  parameter  reduction  is  the  scheme 
which  defines  unloading  and  reloading  elastically.  This  reduces  the  number  of 
parameters  related  to  hysteresis  to  6,  so  the  total  number  of  stored  parameter 
for  the  directional  cracki ng/endochronic  compression  model  is  16. 

RECOMMENDATIONS 

This  investigation  has  shed  considerable  light  on  the  potential  of  the 
endochronic  concrete  model  as  a  practical  analytical  tool.  However,  several 
areas  have  been  uncovered  which  require  investigation.  Consequently,  it  is 
recommended  that  consideration  be  given  to  the  following  areas  for  further 
study: 

1.  There  is  some  discrepancy  between  the  model  and  the  relatively 
meager  data  available  for  loading  paths  parallel  or  nearly  parallel  to  the 


failure  surface.  It  is  recommended  that  the  model  be  tested  against  data 
recently  generated  at  the  University  of  Colorado  which  studies  these  loadings 
i n  some  detai 1 . 

2.  A  problem  seen  in  matching  some  of  the  NMSU  stress  paths  was  the 
excessive  accumulation  of  the  intrinsic  time  increment.  The  excessive  accumu¬ 
lation  was  often  associated  with  stress  paths  which  involved  some  cyclic 
behavior.  It  is  recommended  that  a  detailed  study  be  made  of  the  influence  of 
the  load-unload  criteria  on  the  behavior  of  the  model.  Alternate  criteria 
would  be  a  consideration. 

3.  The  initial  approach  to  strain  rate  effects  made  in  this  study  for 
the  endochronic  model  proves  to  be  quite  promising,  based  on  comparison  with 
the  limited  amount  of  one-dimensional  data.  Comparisons  with  the  limited 
amount  of  multidimensional  data  are  now  desirable,  followed  by  interaction 
between  the  model  and  a  strong  experimental  program. 

4.  The  parameter  reduction  investigation  was  done  without  permanently 
changing  the  endochronic  model  coding.  The  change  should  be  made  permanent 
through  appropriate  changes  in  the  code.  This  would  result  in  a  separate  code 
for  use  when  hysteresis  effects  are  not  needed  in  the  model. 

5.  At  this  point,  consideration  should  be  made  to  expansion  of  the 
model  to  include  the  effects  of  reinforcing  steel.  This  would  make  the  model 
more  efficient  and,  ultimately,  more  accurate  as  more  aspects  of  reinforced 
concrete  behavior  are  introduced.  This  problem  could  be  approached  in  three 
stages.  The  first  stage  involves  a  simple  smearing  of  the  reinforced  steel; 
the  second  would  introduce  aspects  of  steel  concrete  interaction;  and  the 
third  phase  would  be  an  investigation  of  an  anisotropic  model  which  could 
handle  the  effect  of  the  steel  as  efficiently  as  possible. 
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